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ABSTRACT 

The H II region W40 harbors a small group of young, hot stars behind roughly 9 magnitudes 
of visual extinction. We have detected gaseous carbon monoxide (CO) and diatomic carbon 
(C2) in absorption toward the star W 40 IRS la. The 2-0 RO, Rl, and R2 hues of ^^CO at 2.3 
Ilia were measured using the CSHELL on the NASA IR Telescope Facility (with upper limits 
placed on R3, R4, and R5) yielding an Nco of (1.1 ± 0.2) x 10^^ cm~'^. Excitation analysis 
indicates T^in > 7 K. The Phillips system of C2 transitions near 8775 A was measured using 
the Kitt Peak 4-m telescope and echelle spectrometer. Radiative pumping models indicate a 
total C2 column density of (7.0 it 0.4) x 10^^ cm~^, two excitation temperatures (39 and 126 
K), and a total gas density of n ~ 250 cm^'^. The CO ice band at 4.7 jjlvh was not detected, 
placing an upper limit on the CO depletion of (5 < 1%. We postulate that the sightline has 
multiple translucent components and is associated with the W40 molecular cloud. Our data 
for W40 IRS la coupled with other sightlines, shows that the ratio of CO/C2 increases from 
diffuse through translucent environs. Finally, we show that the hydrogen to dust ratio seems 
to remain constant from diffuse to dense environments, while the CO to dust ratio apparently 
does not. 

Subject headings: dust, extinction - infrared: ISM: lines and bands — ISM: abundances — 
ISM: clouds — ISM: individual (W40) — ISM: molecules 



1. INTRODUCTION 

Molecular species are useful diagnostics for under- 
standing the physics and chemistry of the interstel- 
lar medium (ISM), and have been studied in various 
ways for many years. Most molecules are observed in 
dense molecular clouds via rotational emission lines 
in the radio band. These emission studies are in- 
dispensable to our understanding of galactic ecology 
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and the ISM: They yield maps of dense regions, have 
very high spectral (velocity) resolution, and allow us 
to study complex molecules not otherwise observable. 
There are some drawbacks to molecular emission 
studies, however. First, analysis of the line excitation 
is very model-dependent and can lead to significant 
systematic errors. Second, since spatial resolution is 
typically low (and varies from species to species), it 



CO and C2 Absorption Toward W 40 



is sometimes hard to discern whether emission from 
different species comes from the same location within 
a cloud. This makes comparisons among species dif- 
ficult. 

Molecular absorption line observations have some 
important advantages over emission line studies. 
First J absorption from different species observed in 
the same line of sight are more likely to coexist spa- 
tially, thus reducing geometrical ambiguities and al- 
lowing more reliable inter-species comparisons. Sec- 
ond, because the unexcited molecules are observed 
directly, column density determinations are not as 
model-dependent as for emission studies, and hence 
more accurate. Third, some molecular species (no- 
tably H2 and C2) are essentially unobservable in the 
radio due to lack of a dipole moment. These species 
can sometimes be observed via electronic transitions 



in th ; ultraviolet (e.g., Shull & Beckwith 1982 ), or 
through rotational-vibrational transitions in the vis- 



ible or infrared (e.g. Federman et al. 1994; Lacy et 
al. 1J9^). And finally, molecular abundances derived 
from absorption measures can easily be compared to 



line- qf^ight bulk properties (e.g., extinction, polar- 
ization, and atomic abundances) which are also de- 
rived from absorption measures. 

CO and C2 have transitions at 2.3 /xm and 8775 
A, respectively, which can be observed in absorption. 
The primary drawback to studying molecular absorp- 
tion lines in the infrared (IR) is finding a suitable 
background star. The target must be bright in the IR 
yet fortuitously placed behind a substantial amount 
of absorbing material. 

In our effort to find background sources for absorp- 
tion line studies in regions also accessible to emission- 
line studies, we have selected the target W40 IRS 
la, an OB star embedded within the radio source 
W40. This source appears well-suited for our pur- 
poses; Though dim in the visible {V = 15.0) it is 
bright in the IR {K = 5.6), and lies behind about 9 
magnitudes of extinction in V {Ay)- W40 itself is a 
blister-type H II region breaking out of a local molec- 



radio continuum emission, and CO emission. The ex- 
istence of circumstellar dust shells around W40 IRS 
la, 2a and 3a (T ~ 250 - 350 K and M < O.IM©) 
has been suggested based on broad-band IR imaging 
and continuum measurements in the millimeter and 



sub- millimeter ( Smith et al. 1985 ; Vallee &: MacLeod 
199^ ). 

Both CO and C2 are important interstellar 
molecules. CO is the most abundant molecule in 
the ISM after molecular hydrogen, and has a num- 
ber of important roles: The CO emission lines at 2.6 
mm and shorter are not only responsible for cooling 
molecular clouds but also serve as tracers for II2 and 



dense molecular gas (e.g., Magnani &; Onello 1995 
and references therein). C2 is much less abundant 
than CO and has been observed primarily in diffuse 



and translucent clouds (e.g. van Dishoeck & Black 



19891 ; [Lambert, Sheffer, fc Federman 1995| ). Like H2, 



C2 has no dipole moment and hence no pure rota- 
tional spectrum. It is a very useful diagnostic of cloud 
physical conditions such as kinetic temperature, den- 
sity, and radiation field intensity ( |van Dishoeck & 



Black 1982; yan Dishoeck 1984). Both CO and C2 



are important to carbon chemistry everywhere in the 
ISM. 

CO absorption at 2.3 /xm (ti = 2 — 0) has not been 
observed as frequently as the v = 1 - band at 4.7 
/xm, most likely owing to the smaller transition prob- 
abilities. It is, however, easier to observe at 2.3 /xm 
since the thermal background from the Earth's at- 
mosphere is much lower than at 4.7 /xm. Black and 
Willner (1984) and Black et al. (1990) used the lines 
at 2.3 fxm to study the physical conditions and chem- 
istry toward NGC 2024, NGC 2264 and AFGL 2591. 
Recently, Lacy et al. (1994) were able to observe both 
CO and II2 absorption near 2.3 /xm toward NGC 2024 
IRS 2 yielding for the first time a direct measure of 
the H2/CO ratio in a molecular cloud. 

Absorption lines of C2 have been used to address 
a number of problems in the ISM, including: Diffuse 



ular (|;loud, about 400 pc away towards the galactic 



center ( |Zeilik fc Lada 1978| ; [Crutcher fc Chu 1982| ). 
A group of hot stars (most likely B-class) is bathing 
the region with ionizing radiation; W40 IRS 2a (OS 



interstellar cloud chemistry (e.g. Lambert, Sheffer, 
Federman 19951), carbon chemistry in translucent 



clouds (e.g. van Dishoeck fc Black 1989D, molecular 
cloud envelopes (e.g. Federman et al. 1994|), and the 
line of sight structure toward ( Oph (e.g. [Crawford 



2a) s eems to be producing most of the energy (^mith| 1997 ). Federman et al. (1994) provide a good com- 



et al. 1985 ). IRS la is brighter than 2a in the optical 
and near-IR indicating somewhat greater extinction 
toward IRS 2a. Crutcher and Chu (1982) mapped 
i^CO and ^^CO, HCO+, HCN and Ha emission to- 
ward W40, generating a comprehensive kinematic in- 
terpretation of the region. Vallee et al. (1987, 1991, 
1992, and 1994) have done extensive studies on the 
properties of the molecular cloud and its interaction 
with the H II region using C II recombination lines. 



pilation of C2 measurements up to 1992. 

In this paper we report on absorption-line studies 
of CO and C2 toward W40 IRS la. The line of sight 
is discussed in the next section. In Section 3 we de- 
scribe our observations and the results, and in Section 
4 we discuss their implication for the physical state 
of the material in this line of sight. The final section 
contains a brief summary of our conclusions. 
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2. LINE OF SIGHT 

Very little is known about the interstellar sight- 
line to W40 IRS la. Both Crutcher and Chu (1982) 
and Smith et al. (1985) infer Ay ~ 9. The aver- 
age extinction due t o diffuse mate rial over 400 pc is 
about 0.6 mag (see Spitzer 197^ , p. 155). There- 
fore the obscuring material is most likely local to the 
W40 region and more dense than the diffuse ISM. In 
addition, Crutcher and Chu (1982) found a molecu- 
lar ^^CO component in front of the W40 H II region 
with Vlsr = 8 km s'^, Nco - 1-2 x 10^*^ cm'^, and 
log (riT) ~ 4, further suggesting that the sightline 
passes through cold material, most likely associated 
with the neighboring molecular cloud. 

The line of sight intersects at least 3 distinct phys- 
ical regimes. As noted above there appears to be 
some cold foreground material, perhaps associated 
with the nearby molecular cloud. Closer to W40 IRS 
la there must be a photon-dominated region (PDR), 
and closer still, the W40 H II region itself. If W40 
IRS la has a circumstellar dust shell, it is not clear 
what its effect on the CO and C2 spectra might be. A 
warm dust shell might be associated with an elevated 
gas temperature, and hence the appearance of high- J 
molecular rotational lines. 

The H II region is roughly 0.9 pc across and has 
an electron density of 200 cm~'^ (based on results in 
Crutcher and Chu [1982] and a distance of 400 pc). 
The column density of ionized hydrogen associated 
with the H II region should be about 3 x lO^'^ cm~^, 
assuming spherical geometry and that all hydrogen 
is ionized. Using the relation of total hydrogen col- 
umn density to extinction found in Bohlin, Savage, 
&: Drake (1978), the H II region should contribute 
Ay ~ 0.1 mag to the line-of-sight extinction. This, 
of course, assumes a standard gas-to-dust ratio, which 
may not apply to H II regions in general. The con- 
tribution to Ay could increase or decrease depend- 
ing primarily on the nature of any grain-destroying 
shocks which may have passed through the region. In 
either case, the H II region probably does not strongly 
contribute to the total extinction on the line of sight 
to W40 IRS la. In addition, molecules like CO and 
C2 cannot survive in the harsh H II environment, 
so the H II region should not contribute to the line- 
absorption for these molecules either. 

There is almost surely a photon-dominated region 
(PDR) on the line of sight to W40 IRS la, which 
could account for as much as Ay ~ 10, nearly the en- 
tire measured visual extinction (see [Hollenbach 1990] 
for a good discussion). Hydrogen is expected to be 
neutral or molecular throughout the PDR. CO and 
perhaps C2 cannot survive in the regions of the PDR 
closest to the H II region. The effect of dust in PDRs 
is poorly understood, and so a quantitative treatment 



of these regions in general is difficult. 

In summary, CO and C2 absorption lines should 
sample the cold material in the foreground, at least 
part of the PDR, and none of the H II region. A cir- 
cumstellar dust shell, if it exists, would most likely 
only contribute to high- J molecular absorption. The 
visual extinction (Ay ~ 9) should arise almost en- 
tirely in the foreground material (which may be part 
of the local molecular cloud) and the PDR, though 
we note that the PDR could in theory produce all of 
the observed extinction. 

3. DATA AND ANALYSIS FOR ^^CO AND C2 TOWARD 
W40 IRS la 

3.1. C2 Observations at 8775 A 

The observations of W40 IRS la were obtained with 
the Cassegrain echelle spectrograph and the RCA 
CCD camera on the 4-m Mayall telescope at Kitt 
Peak National Observatory on the night of UT 16 
June 1983. Two spectra, each of one hour's dura- 
tion, were recorded in the region of the 2-0 band of 
the Phillips system (A^H^ " X^S+) of C2. The en- 
tire 2-0 band was contained in a single order with its 
center near 8775 A and with a nominal reciprocal dis- 
persion of 0.074 A per pixel at the face of the CCD 
chip. Quartz lamp and thorium-argon spectra were 
obtained for flat-fielding and wavelength calibration. 
An 84 //m wide slit was employed, providing a nomi- 
nal resolution of 0.15 A which in turn corresponds to 
approximately 2 pixels on the chip. 

The two consecutive one-hour frames were co- 
added and averaged. Similarly, two separate flat-field 
frames were co-added and averaged, and the result 
was divided into the averaged stellar frame to produce 
the final photometric spectrum. Both flat-field and 
stellar frames were bias-corrected. The average of 
fifteen bias frames was first subtracted from the flat- 
field and stellar frames, after which a second-order 
bias correction was accomplished by subtracting from 
the average stellar and fiat-fleld frames the mean 
row biases obtained from the masked bias columns 
of each. 

The final spectrum was extracted by collapsing the 
three columns along the slit direction which contained 
the maximum signal. This spectrum is shown in Fig- 
ure 1. Spikes are due to cosmic ray hits and possi- 
bly OH airglow lines. Rotational lines in all three 
branches (P, Q and R) were identified, with the rota- 
tional quantum number J reaching as high as 12 for 
the Q-branch. Unfortunately, there were problems 
with the wavelength calibration at the telescope and 
we were not able to derive a precise measure for the 
radial velocities of the C2 lines. Equivalent widths for 
these lines were determined using standard spectral 
reduction procedures in the NOAO/IRAF package. 
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and column densities for each of the rotational lev- 
els were calculated from the rotational lines using a 
simple Gaussian curve-of-growth (c.f. ppitzer 1978| ). 
Initially, the /-value derived by Erman et al. (1982) 
was used. Recently, the Phillips system transition 
probabilities have been refined (see Lambert, Sheffer^ 
& Fee [erman 1995| and references therein) and the col- 
umn densities have been adjusted to reflect the new 
/-value derived by Lambert, Sheffer, and Federman 
(1995), (1.23±0.16)xl0-3. The data for each line and 
the column densities for each J are shown in Tables 1 
and 2. Several of the lines exhibited effects of satura- 
tion, and an estimate of the Doppler constant (6) was 
derived by requiring that rotational lines originating 
from the same rotational energy level yield the same 
rotational population for that level. This method is 
a simple extension of the doublet ratio method used 
in the analysis of saturated atomic species. The best 
value for the Doppler constant was found to be 1.25 
km s~^. The J = 12 data are not included due to 
high uncertainty: Only the Q12 line was observed 
and it is very weak. 

3.2. ^^CO Observations at 2.3 ^m 

The f = 2 — rotational- vibrational lines for CO 
fall near 2.3 /im, within the K band. Since the os- 
cillator strengths are smaller, the f = 2 — lines are 
not as saturated as their f = 1 — cousins at 4.7 //m. 
In addition, it is somewhat easier to observe at 2.3 
;um, as the sky emission at 4.7 //m is much greater 
and more problematic. Contamination from stellar 
CO absorption should be negligible since W40 IRS 
la appears to be a hot O or B star, which would not 
allow stellar CO to survive ( Crutcher fc Chu 1982 ; 
Smith et al. 198^ ). Circumstellar gas and dust, if it 



exists, could affect the high-J molecular levels. 

Ro-vib transitions at 2.3 /^m for ^^CO were ob- 
served UT 11 June 1994 and UT 21 July 1997 us- 
ing the CSHELL IR spectrometer at NASAs Infrared 
Telescope Facility (IRTF) on Mauna Kea. The IRTF 
is a 3 meter primary, off-axis cassegrain yielding 
f/13.67 at the spectrometer slit. The CSHELL is a 
cryogenically cooled echelle spectrometer with a 256 
X 256 SBRC InSb detector array ( Greene et al. 1993| ). 
We used the 0.5 " slit which gives R ~ 43, 000. Only 
the RO, 1, and 2 transitions were detected (Figure 3) 
with upper limits placed on R3, 4 and 5. A sum- 



maryi of the ubsei vaLiuns is shown in Table 3. — The 
data were obtained and reduced following typical IR 
observing procedures which we briefly summarize be- 
low. 

Dark frames were coadded and subtracted from the 
flat-field image for each wavelength setting. W40 IRS 
la and the standard stars were observed in both the 
"A" and "B" beams (each beam places the spectrum 



in a different spatial location on the detector). Beam 
differences (A-B and B-A) were calculated to elimi- 
nate sky emission, then coadded and normalized us- 
ing the flat field for the appropriate wavelength set- 
ting. Wavelength calibration was achieved using Ar, 
Kr, and Xe lamps with 3 lines per wavelength set- 
ting. One dimensional spectra for W40 IRS la, the 
standard stars, and the calibration lamps were then 
extracted from the detector images using the APALL 
task in IRAF. 

Once the wavelength solution was applied, telluric 
absorption features were identified in both the stan- 
dard stars and W40 IRS la spectra. Telluric lines 
were eliminated from the W40 IRS la spectra using 
the IRAF task TELLURIC, which shifts and scales 
the standard star spectrum before dividing into the 
object spectrum. It is important to note that merely 
dividing the standard star spectrum into the object 
spectrum does not produce the true object spectrum 
without telluric features. To properly remove telluric 
features, one must generate an expected object spec- 
trum, convolve it with the standard star spectrum, 
and compare to the actual observed object spectrum. 
Lacy et al. (1994) give a good discussion of this tech- 
nique. We opted to merely divide, however, as our 
data quality did not warrant more sophisticated tech- 
niques. 

The RO, 1, and 2 lines are shown in Figure 3. These 
lines imply visr = 2 it 2 km s~^, in contrast with the 
8 km s~^ foreground molecular material seen in emis- 
sion ( Crutcher fc Chu 1982| ). Since emission line data 
sample a large, beam-averaged area, and absorption 
lines just a pencil beam, we do not necessarily expect 
velocities derived from both to agree. We merely note 
that they are not wildly different. 

R = 43, 000 corresponds to a resolution of Aii ~ 7 
km s-^ The RO, 1, and 2 hues all have FWHM ~ 8 
km s~^ and hence are not well-resolved. After contin- 
uum normalization, each line was directly integrated 
for equivalent width. Errors reflect continuum place- 
ment ambiguity and are of high confidence {2a). All 
the widths (along with transition information and er- 
rors) are given in Table 4. Widths for the R3, R4 and 
R5 lines are 2a upper limits based on the noise in the 
continuum at the expected line position. 

The equivalent width of any absorption line is de- 
pendent on the column density of the species and, 
if saturated, the velocity parameter b (c.f. ^pitzei) 
1978|) . We generated a model curve of growth (COG) 
using b = 1.25 km s~^, as determined from the C2 
lines. Each ^^CO equivalent width was fit to the 
COG independently and the column density for each 
rotational level {Nj) is given in Table 5. The R2 
through R5 lines were clearly optically thin while the 
RO and Rl lines were more optically thick, falling 
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the transition from the hnear portion to the satu- tor for the radiation field in the far-red (van Dishoeck 



on 

rated part of the COG. If the 6- value for CO is greater 
than 1.25 km s~^, then the RO and Rl lines become 
more optically thin and their abundances drop by 2 
and 0.5 xlO^ cm~^ respectively. Summing the abun- 
dances for each line gives a total column density of 
Nco = (1.1 ±0.2) X 10^^ cm~^ (assuming ^^CO to be 



fc Black 1982) 



Total molecular abundance and excitation tem- 
peratures were calculated from models which fit a 
Maxwell-Boltzmann population distribution to the 
observed rotational level abundances. As had been 



found for other relatively dense clouds (Lutz & 



the d jminant isotope). Depending on the excitation Crutcher 1983 ), the rotational abundances could not 



of the higher J levels, this value could be slightly too 
small, and if 6 > 1.25 km s~^, then it would be too 



high. Our column density for CO is nearly identical 



to that inferred by Crutcher &; Chu (1982) for the 
foreground, Nco ~ 1-2 x 10^ 



ll8 pjj^-2 



3.3. CO Ice Band at 4-7 fnn 

In an effort to detect the CO ice band at 4.7 /um we 
made observations with moderate spectral resolution 
at the United Kingdom IR Telescope (UKIRT) using 
the CGS4 on 19 August 1998. All observations were 
carried out while nodding along the slit to remove 
atmospheric emission. The total integration time on 
W40 IRS la was 11.2 minutes at an average airmass 
of 1.45. The spectrum for W40 IRS la was ratioed by 
BS 7236 (B9V) to remove telluric absorption features. 
The S/N is ~ 44 at 4.67 /xm and no CO ice feature is 
apparent. An upper limit for the optical depth of the 
band is T 4.67 < 0.02 {2a), implying N cojlce) < lO^*' 



cm ( Tielens fc Allamandola 1987 ). The depletion 



of CO into icy mantles for the W40 IRS la sightline 
must be less than 1 %. 



4. DISCUSSION 



4.1. Excitation Conditions and the Foreground 
Cloud Physical Properties 



Welhave three basic diagnostics for the temperature 



and density toward W40 IRS la. The ^^CO emis- 
sion data imply log (nT) ~ 4, much low er than the 
nearby molecular cloud, log (nT) ~ 6.5 ( prutcher & 
Chu 1982). Excitation conditions can also be derived 
from the rotational populations of CO and 02- 

The rotational level populations of the C2 Phillips 
system {v = 0, J) are attained via radiative pump- 
ing. Lifetimes of these levels are so long that collisions 
and upward electronic transitions are the most impor- 
tant depopulation mechanisms. Hence the rotational 
populations reflect the competition between pumping 
and collisions and are non-thermal in general ( |van| 
Dishoeck fc Black 1982 ). For densities greater than 



be fit with a single temperature distribution, presum- 



ably as a result of radiative pumping ( Chaffee et al 
198C; van Dishoeck fc Black 1982| ). Consequently, 



we fit these data with a two-temperature model used 
by Lutz and Crutcher (1983): J-levels 4 through 10 
were fit with an excitation temperature {Te^) of 126 
K, which we associate with the effects of radiative 
pumping. After correcting the populations in J = 
and 2 for the contributions from the 126 K pop- 
ulation distribution, we derived a J=2/J=0 excita- 
tion temperature of 39 K, which we associate with 
the thermal distribution of the gas. The resulting 
total column density of C2 towards W40 IRS la is 
(7.0 ± 0.4) X 10^^ cm~^, based on a Boltzmann dis- 
tribution for both temperatures. This is the highest 
column density of C2 yet seen in absorption. Figure 
2 shows the final fits. 

In comparing the results for W40 IRS la to the ra- 
diative pumping models calculated by van Dishoeck 
(1984), we found that they are best characterized by 
her model with a kinetic temperature of 40 K and a 
radiation parameter, ^^^ of 5.27 x 10~^^. Assuming 

Ir = 1, cr ~ 2 X 10"^*^ cm2 ( |van Dishoeck fc Black 
1989| ), and that hydrogen (H I + H2) is the only im- 
portant collision partner, we get Ufj ~ 250 cm~^ for 
the W40 IRS la sightline. Since the value of a is 



not known to better than a factor of 2 (van Dishoeck 
Black 1989| ), our value for uh is equally impre^ 



cise. In addition, for an enhanced (or depleted) radi- 
ation field, the estimated collision (hydrogen) density 
would scale accordingly. 

Rotational transitions of CO are primarily excited 
by collisions with hydrogen. De-excitation can occur 
via collisions with hydrogen or by spontaneous line 
emission. The critical density at which collisions be- 
gin to overtake spontaneous emission is around 3000 
cm^^. More detailed studies of the excitation, pho- 
todissociation, and chemistry have been conducted 
by van Dishoeck & Black (1988) and Warin, Benay- 
oun, & Viala (1996). In general it is found that the 
rotational populations of CO are sub-thermal except 



100 c m"'^ the J = a.nd J = 2 levels are verv nearlv ^^ ^^^ ^^^^ ^^^ lev els in dense cases (Warin, Benay 



thpr JiRl OthprwisP thP pnpiilRtinns dppPTid on thp oun, fc Viala 199^ ). Hence the assumption of LTE 



thermal temperature, T, and the radiation param- 
eter, Y^) where ric is the collision partner density 
(usually n(II) + n(H2)), a is the effective cross section 
for collisional de-excitation, and Ir is the scaling fac- 



will almost always produce excitation temperatures 
which underestimate the actual thermal temperature 
{Tex < Tkin)- Recent work by Wannier, Penprase, & 
Andersson (1997) suggests that the dominant form of 
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CO excitation in diffuse and translucent clouds can 
be line emission from nearby molecular clouds, if the 
clouds have similar velocity vectors. 

As a start, we constructed a Boltzmann plot for 
the ^^CO data (Figure 4). Fitting the temperature 
to all lines gives T^x = 7 K. This value is similar 
to that found by Cruther & Chu (1982), but it is 
important to note that the excitation temperatures 
for ^^CO and ^^CO are different in general ( War in J 
Bena youn, fc Viala 1996| ). Since the thermal tem- 
perature is greater than the CO rotational tempera- 
ture in general, T^in > 7 K, which is consistent with 
the C2 analysis. Warin, Benayoun, & Viala (1996) 
have constructed models of CO excitation for the 
dense, translucent, and diffuse cloud regimes, includ- 
ing UV photodissociation. These models show very 
clearly that the excitation temperature is much lower 
than the thermal temperature for diffuse and translu- 



cent clouds. The excitation of CO in dense clouds is 



nearly thermal and T^x = Tf^m- If we assume that the 
cloud(s) toward W40 IRS la are dense, then T^in — 7 
K. For diffuse and translucent clouds, the level pop- 
ulations observed can be scaled (see Section 4.2 of 
Warin, Benayoun, & Viala [1996]) to derive "pseudo- 
LTE" rotational populations, i.e. representative LTE 
populations for the actual thermal temperature of the 
gas. A slope can be fitted to these populations and a 
thermal temperature inferred. If the cloud(s) toward 
W40 IRS la are translucent, then the RO population 
is reduced by ~ 1.75 (changes in Rl and R2 happen 
to be negligible). The resulting populations imply 
Tfcjn ~ 9 K (see Figure 4), which still does not agree 
with the C2 analysis. It is not clear, however, that 
the models constructed by Warin, Benayoun, & Viala 
(1996) apply to our line of sight: The UV field is prob- 
ably higher than normal near W40, and it seems very 
likely that there are multiple clouds on the sightline. 
Since the cloud(s) on the W40 IRS la line of sight 
are very near the W40 molecular cloud, and the LSR 
velocities of both are similar (5 and 8 km s~^), the ex- 
citation for ^^CO may be at least partially radiative. 
Comparing to models constructed by Wannier, Pen- 
prase, &: Andersson (1997) it is apparent that our 
Tex = 7 K is degenerate: It can be accounted for 
by many combinations of radiative excitation, colli- 
sion partner density, and kinetic temperature. The 
C2 excitation conditions can help constrain those for 
^^CO. If we assume T(fcj„ ~ 40 K, the density of the 
gas can range from to 450 cm~^, depending on the 
efficiency of radiative excitation. Collisional excita- 
tion becomes dominant as n approaches 450 cm^^ 
(which we use as an upper limit). If we further as- 
sume n ~ 250 cm~^ (as indicated by the C2 data), 
then the excitation of CO can only be explained by 
collisional and radiative processes combined. 



A summary of the W40 IRS la sightline cloud phys- 
ical properties based on the CO and C2 data in this 
paper as well as the CO emission study by Crutcher 
and Chu (1982) is given in Table 6. As discussed in 
Section 2, we assume that nearly all of the absorp- 
tion on this line of sight arises in one cloud complex 
(single or multiple components) local to the W40 re- 
gion, since it is only 400 pc distant (i.e., the cloud(s) 
cannot be diffuse). The physical conditions are best 
determined by C2: ra ~ 250 cm~^ and T ~ 40 K. The 
CO emission and absorption values and limits agree 
very well. The cloud(s) are clearly not dense enough 
to be considered molecular, but could be considered 
translucent. Translucent lines of sight typically have 

Ay = 2 5 and can be studied via both absorption 

lines (UV, optical, and/or IR), and radio emission 
lines. In addition, translucent material is expected in 



the outer envelopes of molecular clouds (van Dishoeck 
fc Black 1989|) . We postulate that the W40 IRS la 



sightline is composed of multiple translucent compo- 
nents associated with the W40 molecular cloud. 

4.2. Molecular Abundances 

CO and C2 abundances have been determined 
jointly for a number of sightlines, a sample of which 
is shown in Table 7. W40 IRS la is one of few sight- 
lines allowing a direct comparison of CO and C2 in 
absorption. The abundance ratio we derive is CO/C2 
= 1600 ± 600 {2a). The depletion of CO onto dust 
is negligible (< 1%), in view of our failure to detect 
CO ice, but the C2 depletion has not been assessed. 
Note that the CO/C2 ratio increases from diffuse to 
translucent and molecular regimes indicating (to first 
order) that the formation/destruction rates favor CO 
over C2 as cloud type changes. This may merely be 
due to self-shielding of CO, but other processes may 
also be in action: The data are not yet precise enough 
to tell. 

Assuming all hydrogen is molecular, we can esti- 
mate the amount of H2 on the W40 IRS la line of 
sight from Nco- Using an H2/CO ratio of 3700l:^7[![J 
(based on the direct comparison of H2 and CO IR ab- 
sorption lines toward NGC 2024 IRS 2 {Ay = 21.5 ± 



5), ILacy et al. T994|) , we find Nh^ = 4.4+^;^ x 10' 



1+3.8 



21 



cm~ . The H2:C0 ratio is derived from a direct com- 
parison of weak absorption lines and hence should be 
reliable, but is based on only one data point: NGC 
2024 IRS 2. If this sightline is abnormal in any way, 
then the ratio does not necessarily apply to other lines 
of sight such as W40 IRS la. 

4.3. Dust Indicators 

The Ay calculated by Crutcher and Chu (1982) and 
Smith et al. (1985) assumes a "normal" extinction 
law {Rv ~ 3), which may be incorrect. The ratio 
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of visual to selective extinction, Ry, is a grain size 
distribution indicator ( [Cardelli, Clayton, &: Mathis 
1989|]|: Ry < 3 implies an abundance of small grains 
compared to the typical size distribution, whereas 
Ry > 3 implies that larger grains dominate the 
extinction. As interstellar clouds collapse the dust 
grains tend to agglomerate, eliminating the smaller 
particles ( Jura 1980| ). If the line of sight toward W40 
IRS la does indeed sample the edge of a molecular 
cloud, then we would expect a population of larger 
grains and 3 < Ry < 5. Smith et al. (1985) found 
{B — V) = 2.2 for W40 IRS la and assuming it is an 
OB star, Es-y — 2.5. For Ry = 3 — 5, we get a range 
in Ay of 7.5 to 12.5, which agrees with the previous 
work by Crutcher and Chu (1982) and Smith et al. 
(1985). 

Our calculated value of Ay can be used to pre- 
dict the amount of CO expected toward W40 IRS la. 



As dscussed in Section 2, however, the extinction 



and Federman (1995)), about half of the abundance 
on the W40 IRS la sightline. Recently, McCall et al. 
(1998) measured CO IR absorption toward Cyg 0B2 
No. 12 yielding Nco = 2 x 10^^ cm'^, a factor of 60 
less than what we find for W40 IRS la over nearly 
the same total visual extinction. The 4.7 fixa CO ice 
feature is not seen on this line of sight; therefore, de- 
pletion of CO onto ice mantles cannot readily explain 
the discrepancy. This clearly shows that the gaseous 
CO to dust ratio changes from diffuse to denser en- 
vironments. 

As an interesting side note, we have compared the 
total hydrogen column density (assuming all of it to 
be molecular) for both W40 IRS la and NGC 2024 
IRS 2 to the hydrogen/reddening correlations found 
by Bohlin, Savage, and Drake (1978) and Dickman 
(1978) (Figure 5). The Ee-y for NGC 2024 IRS 
2 assumes Ay = 21 ( |Lacy et al. 1994 ; piang, PeF 
rier, fc Lena 198^ and Ry ranging from 3 to 5. The 



and CO absorption may not be well-correlated in the 
PDR, leading us to slightly overestimate Nco based 
on Ay. Conversely, most studies of the CO/ Ay ra- 
tio have not addressed the depletion of CO into ice 
mantles (which may be as high as 40 %, [Chiar et"^ 



1995|j. Since CO is almost entirely in the gas-phase 
for the W40 IRS la line of sight, these studies will 
underestimate the predicted CO column density. Us- 
ing radio maps, spectroscopic data, and star counts 
in the Taurus and p Oph clouds, Frerking, Langer, 
and Wilson (1982) found: 



NiC'^'O) = 1.7 X 10^^(^v' - 1-3) 



(1) 



for N in cm ^ and 4 < Ay < 21. This rela- 
tion predicts Nco = (0.5 - 1.0) x 10^^ cm'^ (with 



N{CO)/N{C'^^0) = 490) for our line of sight, very 
nearly the same as we have measured. Interestingly, 
it may be that the PDR effect nearly balances the 
depletion eff'ect. 

There is no strong correlation between C2 and 
Eb-v ( van Dishoeck fc Black 1989| ) though other 
dust indicators have not been investigated. The lack 
of correlation may be due to the fact that many of the 
stars included in C2 absorption studies to date have 
been distant supergiants where some of the sight- 
line extinction is caused by diff'use clouds, which lack 
abundant 02- 

Comparing the W40 IRS la sightline to others in 
Table 7, it is apparent that CO and C2 generally in- 
crease with Ay but the data are too scattered to 
draw any strong conclusions. 

The measurement of CO and C2 absorption toward 
Cyg 0B2 No. 12, the classic diffuse cloud line of sight 
with Ay ~ 10, is quite interesting. Lutz and Crutcher 
(1983) found Nc^ = (3.0 ±0.2) x 10^"^ cm'^ (adjusted 
for new Phillips system /-values in Lambert, Sheffer, 



correlations ( "intercloud" and "cloud") from Bohlin, 
Savage, and Drake (1978) are based on La absorp- 
tion for lightly reddened sightlines with Es-y < 0.6, 
while the relation from Dickman (1978), derived from 
CO emission, is good to Es-y ~ 3. In Figure 5 we 
have extrapolated out to Es-y = 10. The data from 
W40 IRS la and NGC 2024 IRS2 agree with the ex- 
trapolated relations surprisingly well, within factors 
of 2 or so. 

5. CONCLUSIONS 

We have used ^^CO and C2 IR and visible absorp- 
tion lines to investigate the line of sight toward W40 
IRS la. The C2 data were obtained at the 4 m May- 
all Telescope at KPNO and the ^^CO data from the 
CSHELL on the IRTF. The CO ice band at 4.7 ^m 
was not detected. This sightline is clearly much more 
dense than diffuse, based on Ay ~ 9 and a distance 
of only 400 pc (we calculate a range in Ay of 7.5 
to 12.5, assuming a population of large grains). Our 
^^CO and C2 data show: 

1. The C2 excitation conditions, T ~ 40 K and 
n ~ 250 cm~^, agree with limits determined from 
CO emission and absorption. Since the Ay is large 
and both absorption and emission measures are avail- 
able, we postulate that the W40 IRS la sightline has 
multiple translucent components. 

2. The non-detection of CO ice indicates a CO de- 
pletion oi 5 < 1%. 

3. Comparing to other sightlines in Table 7, we find 
an overall increase in Nco and Nc2 with increasing 
Ay. The data, however, are too scattered to draw any 
further conclusions at this point. The ratio of CO to 
C2 appears to increase from diffuse to translucent and 
molecular sightlines probably due to the self-shielding 
of CO. 
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4. The column density of CO toward W40 IRS la is 
~ 60 times that found for Cyg 0B2 No. 12, the clas- 
sic diffuse line of sight ( McCall et al. 199§| ), despite 
very similar Ay values. This is not a depletion effect, 
and suggests that the CO-to-dust ratio changes from 
diffuse to dense environments. 

5. Finally, the relationship of hydrogen column 
density to interstellar reddening ( Bohlin, Savage,^q 
Drak{! 1978; Dickman 1978| ) was found to be roughly 
consistent with recent data out to Eb-v ~ 10 (Figure 
5). Contrary to the CO-to-dust ratio, the hydrogen- 
to-dust ratio appears to be valid from diffuse to dense 
regimes. 

This research was supported by NASA Graduate 



Student Research Program grant (NGT5-50032) to 
R. Y. Shuping and NASA grant NAG5-4184 to T. P. 
Snow. We would like to thank the staff and opera- 
tors at KPNO, the IRTF (J. Rayner), and the UKIRT 
(T. Kerr) for their help and tireless duty through the 
night. Many thanks to J. Black and the referee S. 
Federman, whose comments and suggestions greatly 
improved the content and quality of this paper. B. 
L. Lutz would like to thank K. Sheth who carried out 
the curve-of-growth analysis and temperature fits for 
the C2 lines. R. Y. Shuping would like to thank B-G 
Andersson, J. H. Lacy, and D. Jansen for their helpful 
input, and also D. E. Schutz for useful conversations 
and support. 



REFERENCES 



Black, J. H. & WiUner, S. R 1984, ApJ, 279, 673 
Black, J. H., van Dishoeck, E. F., Willner, S. P., 



Woods, R. 
Drake, J. F. 1978, ApJ, 224, 



C. i990, ApJ, 358, 459 
Bohlin, R. C, Savage, B. D 

132 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 

245 
Chackerian, C. & Tipping, R. H. 1983, J. Molec. Spectrosc, 

99 431 
Chafl'ee, F. H., Jr., Lutz, B. L., Black, J. H., Vanden Bout, P. 

A., & Snell, R. L. 1980, ApJ, 236, 474 
Chiar, J. E., Adamson, A. J., Kerr, T. H., & Whittet, D. C. B. 

1995, ApJ, 455, 234 
Crawford, I. A. 1997, MNRAS, 290, 41 
Crutcher, R. M. & Chu Y. 1982, Regions of Recent Star 

Formation, 53, Roger, R. S. & Dewdney, P. E. eds. 
Dickman R. L. 1978, ApJS, 37, 407 
Erman, P., Lambert, D. L., Larsson, M., & Mannfors, B. 

1982,ApJ 253, 983 
Federman, S. R., Strom, C. J., Lambert, D. L., Cardelli, J. A., 

Smith, V. v., & Joseph, C. L. 1994, ApJ, 424, 772 
Federman, S. R., Lambert, D. L., van Dishoeck, E. F., 

Andersson, B-G, & Zsargo, J. 1999, in prep 
Frerking, M. A., Langer, W. D., & Wilson, R. W. 1982, ApJ, 

262 590 
Grede'l, R., van Dishoeck, E. F., de Vries, C. P., & Black, J. H. 

1992, A&A, 257, 245 

Gredel, R. , van Dishoeck, E. F., & Black, J. H. 1994, A&A, 

285, 300 
Greene, T. P., Tokunaga, A. T., Toomey, D. W., & Carr, J. S. 

1993, Proc. SPIE, 1946, 313 

HoUenbach, D. J. 1990, ASP Conference Series Vol. 12, The 
Evolution of the Interstellar Medium, L. Blitz ed., p. 167 



Jiang, D. R., Perrier, C, & Lena, P. 1984, A&A, 135, 249 

Jura, M. 1980, ApJ, 235, 63 

Lacy, J. H., Knacke, R., Geballe, T. R., & Tokunaga, A. T. 

1994, ApJ, 428, L69 
Lambert, D.'l., s'heffer, Y., & Federman, S. R. 1995, ApJ, 438, 

740 
Lutz, B. L. & Crutcher, R. M. 1983, ApJ, 271, LlOl 
Magnani, L. & Onello, J. S. 1995, ApJ, 443, 169 
Mantz, A. W. & Maillard, J.-P. 1974, J. Molec. Spectrosc, 53, 

466 
McCall, B. J., Geballe, T. R., Hinkle, K. H., & Oka, T. 1998, 

Science, 279, 1910 
ShuU, J. M. & Beckwith, S. 1982, ARA&A, 20, 163 
Smith, J., Bentley, A., Castelaz, M., Gehrz, R. D., Grasdalen, 

G. L., k Hackwell, J. A. 1985 ApJ, 291, 571 
Spitzer, L.Jr. 1978, Physical Processes in the Interstellar 

Medium, (New York: Wiley) 
Tielens, A. G. G. M. & AUamandola, L. J. 1987, Physical 

Processes in Interstellar Clouds, eds. G. E. Morfill & M. 

Scholer^ p. 333, Dordrecht: Reidel 
Vallee, J.P. 1987, A&A, 178,237 

Vallee, J. P. & MacLeod, J. M. 1991, A&A, 250, 143 
Vallee, J. P., Guilloteau, S. & MacLeod, J. M. 1992, A&A, 266, 

520 
Vallee, J. P. & MacLeod, J. M. 1994, AJ, 108, 998 
van Dishoeck, E. F. 1984, PhD Thesis, Leiden University 
van Dishoeck, E. F. & Black, J. H. 1982, ApJ, 258, 533 
van Dishoeck, E. F. & Black, J. H. 1988, ApJ, 334, 771 
van Dishoeck, E. F. & Blacky J. H. 1989, ApJ, 340, 273 
Wannier, P., Penprase, B. E., & Andersson, B-G 1997, ApJ, 

487, L165 
Warin, S., Benayoun, J. J., & Viala, Y.P. 1996, A&A, 308, 535 
Zeilik, M. & Lada, C. J. 1978, ApJ, 222, 896 



Shuping et al. 



Table 1 
Observed C2 Phillips System Lines 



Transition 



Rest Wavelength (A) 



Equivalent Width (mA) 



RO 
R2 
Q2 
P2 
R4 
Q4 

P4 + q 

R6 
Q6 
P6 
P8 

QIO 
Q12 



8757.7 
8753.9 
8761.2 
8766.0 
8751.6 
8763.7 
8773.3 
8750.8 
8767.7 
8782.3 
8792.6 
8780.1 
8788.5 



37.0 ±3.0 
52.5 ±1.5'" 
55.0 ± 5.0 
18.0 ± 3.0 
42.5 ± 7.5^ 
36.0 ±3.0 
43.0 ± 4.0 
22.8 ±3.0 
31.5 ±4.5 
12.5 ±2.5 
9.5 ±2.0 
17.5 ±2.5 
6.5 ±2.0 



''R2 and RIO are blended 
R4 and R8 are blended 



Table 2 
C2 Column Densities 



J 



Ej/k (K) 



Nj (xlO 



13 



cm 




2 

4 



10 



0.000 
15.635 

52.114 
109.430 
187.572 
286.526 



5.4 ±0.8 
19.8 ±2.0 
11.2 ±1.6 
9.2 ±0.7 
5.6 ±0.7 
4.1 ±0.4 



TotaP 



70 ±4 



^Based on Boltzmann Analysis for both 39 and 126 K 
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Table 3 
CSHELL ^^CO Observations Summary 



Wavelength Range (fim) ^^CO Lines Int. Time (min) S/N" 



Standard Stars 



2.334-2.340 
2.340-2.346 



R3, R4 and R5 



25 
25 



90 BS 6714, 7110 and 5 Aql 
90 Spica 



^S/N for the regions surrounding the applicable lines. 



Table 4 
Observed ^^CO Lines {v 



0) 



Transition 


Rest Wavelength 


i^m)^ 


/-value (xlO-8)^ 


Equiv 


. Width (mA) 


RO 


2.34530523 




8.78 




130 ± 20 


Rl 


2.34326929 




5.89 




100 ± 10 


R2 


2.34127497 




5.33 




60 ±17 


R3 


2.33932336 




5.11 




<20 


R4 


2.33741273 




5.00 




<20 


R5 


2.33554435 




4.94 




<20 














^IMantz & Maillard 1974 










'^hackerian and Tipping 1983| 





Table 5 
^^CO Column Densities 



J 


Ej/k (K) 


Nj (xlO^^ cm-2) 








5.0 ±2.0 


1 


5.532 


4.0 ±1.5 


2 


16.60 


2.2 ±0.7 


3 


33.19 


<0.8 


4 


55.32 


<0.8 


5 


82.98 


<0.8 


Total 




11±2 
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Table 6 
W40 IRS lA Line of Sight Physical Properties 



Species 


Vlsr (km s ^) 


b (km s-i) 


iV (cm-2) 


Te. (K) 


nH%cm-^) 


log {nT) 


C2 

"CO Em.'' 

i^CO 


8 

2±2 


1.25 


7.0 ±0.4 X 10^'' 
1.2 X 10^® 

1.1 ±0.2 X 10^* 


39 and 126'' 

> 7 


-250 
<450'' 


^4 



"Crutcher & Chu 1982 



T — 126 K component due to radiative pumping. 
"H±H2 
''Assuming T ~ 40 K, see Section 4.1 
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Table 7 
SiGHTLiNES With CO and C2 Observations 



Sightline 


Nco (cm "^Y 


iVc2 (cm-2)'' 


CO/C2 


Cloud Type" 


Av 


References'* 


W40 IRS la 


1.1 ±0.2(18) 


7.0 ±0.4(14) 


1600 ± 600 


T 


7.5 - 12.5 


This Work 


Cyg OB2 #12 


2(16) 


3.0 ±0.2(14) 


67 ±5 


D 


-10 


1, 2 


HD94413 


0.7-1.2(16)" 


2.8(13) 


250-430 


T 


2.4 


3 


HD 154368 


0.6-1.5(16)" 


4.6(13) 


130-330 


T 


2.5 


3 


HD 169454 


0.55-1.8(16)" 


5.6(13) 


100-320 


T 


3.3 


3 


Per 


1.1(15) 


1.8(13) 


60 


D 


0.9 


4 


HD27778 


2.5(16) 


3.0(13) 


830 


T 


1.2 


4 


p Oph A 


1.9(15) 


2.1(13) 


90 


T 


1.4 


4,7 


COph 


2.3(15) 


1.79 ±0.06(13) 


130 


D/T 


0.96 


4 


20 Aql 


3(15) 


4.2(13) 


71 


D/T 


0.99 


4 


HD207198 


2.6(15) 


2.4(13) 


110 


T 


1.9 


4 


HD21483 


1.0(18) 


7.4(13) 


100 - 10^ 


T 


1.7 


4 


C Per 


1.2(15) 


2.8(13) 


43 


D/T 


1.0 


4 


XPer 


5.0(15) 


4.2(13) 


120 


T 


1.4 


4 


HD 26571 


6.0(16) 


8.8(13) 


680 


D/T 


0.9 


4 


AE Aur 


1.3(15) 


4.6(13) 


28 


T 


1.6 


4 


HD 110432 


1.0(15) 


2.4(13) 


42 


T 


1.2 


4 


TT ScO 


1.0(12) 


< 1.0(12) 


> 1 


D 


0.19 


5 


P^ Sco 


1.2(13) 


< 1.0(12) 


> 12 


D 


0.62 


5 


LU^ Sco 


4.0(13) 


< 2.0(12) 


> 23 


D 


0.68 


5 


xOph 


3.0(14) 


2.8(13) 


11 


T 


1.2 


4 


9 Cep 


1.7(13) 


1.0(13) 


2 


Df 


1.5 


4 


HD 210121 


3.0(15)" 


5.2(13) 


58 


T (High Lat.) 


~ 1 


6 


A Cep 


1.4(15) 


1.4(13) 


100 


Df 


1.8 


4 



Note.— X(y) = X x 10^ 

^AU column densities derived from UV or IR absorption lines except where noted. 

''Where appropriate, the column density for C2 has been adjusted to account for the new /-value derived 
by Lambert, Sheffer, and Federman (1995). 

"Cloud Types: D - Diffuse, T - Translucent 

** References: 1 - JMcCall et al. 1998| 2 - JLutz fc Crutcher 198^. 3 - [Credel. van Dishocck. fc Black 1994 
4 - Federman e t al. 1 994 | and references therein, 5 - [Lambert, Sheffer, fc Federman 1995| , 6 Grcdcl et al. 



1992| , 7 - [Federman et al. 1999 



"From mm emission lines 
Possibly more than one cloud. 
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Fig. 1.— C2 Phillips system toward W40 IRS la. 
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Figure 2 
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Fig. 2. — Boltzmann plot for rotational levels of C2: Two-temperature model, 39 K (solid) and 126 K (dashed). 
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Figure 3 
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Fig. 3. — CO lines at 2.3 ^m toward W40 IRS la. The continuum for Rl and R2 is offset for clarity. 
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Figure 4 
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Fig. 4. — Boltzmann plot for rotational levels of CO. Squares are observed column densities from Table 5 with Tex = 7 K 
(solid line). Triangles are pseudo- LTE populations assuming material on the line of sight is translucent in nature with temperature 
fit {Tex ~ 9 K, dashed line). See section 4.1 for discussion. 
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Figure 5 
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Fig. 5. — Correlation of reddening and total hydrogen column density. Solid line is "intercloud" relationship extrapolation from 
Bohlin, Savage, and Drake (1978); dashed line is the "cloud" extrapolation; dotted line is relationship found by Dickman (1978). 
The total H column density for both W40 IRS la and NGC 2024 IRS 2 assumes all hydrogen is molecular. See Section 4.3 for 
further discussion. 



